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The operator method is used to construct the solutions of the problem of the polaron in the 
strong coupling limit and of the helium atom on the basis of the Hartree-Fock equation. E — 
—0.1085128052a 2 is obtained for the polaron ground-state energy. Energies for 2s- and 3s-states 
are also calculated. The other excited states are briefly discussed. 



I. INTRODUCTION 



The operator method (OM) of the approximate solution of the Schrodinger equation was suggested in the works of 
Feranchuk and Komarov [1-3]. In agreement with the results obtained in these and subsequent works (see, for example, 
h-H ' refs. [4-12]), the solution in the zeroth approximation of the OM gives quite a simple and universal algorithm 
for obtaining an approximation, which is uniformly applicable in a wide range of variation of parameters of the 
t-H 1 Hamiltonian. Another advantage of the OM is connected with the possibility of regularly calculating the correction 
to the zeroth approximation. Here, in any order the needed calculations are reduced to a simple algebraic procedure 
of expressing the product of excitation creation and annihilation operators in the normal form, and this essentially 
simplified the use of computers. 
1 As shown by the results of many works the OM is very effective for solving various physical problems. However, the 
C*") , question of the application of this method in solving non-linear equations, which are often met in atomic and solid 
O^l ■ state physics, remains unclear at present. 

In this work we investigate in two examples the possibility of using the OM for solving non-linear equations such 
as that of Hartree-Fock: 

(i) in an example of the problem of a polaron in the regime of strong electron- phonon coupling (section II) ; 
' (ii) in an example of the problem of a helium atom (section III). 

It should be noted that such calculations have been carried out on the basis of various methods. But here, we want 
to obtain the exact solution by a numerical method. Comparison of approximate solutions with the exact numerical 
solution enables us to estimate not only the accuracy of the approximate solutions, but also the usefulness of their 
g further study. 

o ■ 

^ 1 II. THE OPERATOR METHOD IN THE POLARON PROBLEM 

. • n 

The problem of a polaron in the strong coupling regime has become of interest since Emin ||13| proposed the 
mechanism of bipolaron Bose-Einstein condensation to explain high-T c superconductivity, and on the other hand, 
the bipolaron is known to exist in this strong coupling regime E3. The Hamiltonian of a non-relativistic particle 
(electron) interacting with the quantized scalar field of lattice oscillations in an ionic crystal has the form (see, for 
example, refs. [6,7,15,16] and reviews [17,18]): 

6 = + E c " + 4f E ( A ^ r ^ + Ke-^cA . (i) 

Here, and Ck are respectively the operators of creation and annihilation of phonons with frequency Wk and mo- 
mentum k, Ak are the Fourier components of the source density and f2 is the volume of the system. In the theory of 
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the so-called optical polaron, the frequency of phonons is assumed to be independent of their wavevector. Usually, it 
is assumed that 

iuj {Ana) 1 / 2 

^k = ^0, gA k = — — 7TT7, 

k (2mwo) ' 

where a is a dimensionless electron-phonon coupling constant. 

For the application of the OM it is conducive to cross to the four-dimensional oscillator representation using the 
transformation [|l9| 

x\ = C(°"A) s iCt 
X = arg(£i), 

where a\ (A = 1,2,3) are Pauli matrixes, £ s (s = 1,2) are regarded as spinor components. The physical meaning of 
"extra" variable \ nas been treated in many works (see, for example, refs. [20,21] and references therein). 
Let us define operators 

= \\% it* + ~ J , b a (u) = <l% ( C ' 



2 v wflft/ V 2 V udi 

Here, the positive parameter cj will be defined later. Operators (2) satisfy the commutation relations 

[a 3 (w),af(u)] = S st , [b s (uj),bf(uj)} = S st 

(we have written only non-zero commutators). The possibility of using the algebraic method is conditioned by the 
fact that, all "physical" operators can be expressed through the following 15 operators 

M = a s b s , M+=a+b+, N = afa s +bfb s , (3) 

m x = {<7\) st a t b Sl mX = (cr A ) st a+&^~, n a x = {t7\) st a^a t , n\ = (a\) st bfb Sl 

which form a closed algebra so(4, 2) (see ref. |l^]). Here and henceforth, we omit for brevity the parameter uj in 
expressions of the operators. For further use we rewrite some operators into this representation as follows 

xx = ^(m+ + m x +n a x + n\), r = = i-(M+ + M + N + 2), 

f) 2 LlI 

rA = — — — — — = —(2 + N — M — M+), 

where the operators TV, M, M + form the subalgebra so(2, 1) and satisfy the following relations 

[M, M + ] = 2 + N, [M, 2 + N] = 2M, [2 + N, M+] = 2AI+ . (4) 

In the four-dimensional oscillator representation (£-space) the equation for eigenvalue problem H |\&) = E\^) has the 
form H| 

L\V) =r(H-E)\^) = 0, (5) 

which, in terms of the operators (3), rather simplifies the algebraic calculation of the matrix elements. It is easy 
to see that this representation is especially useful for the bound polaron problem and leads to the equation without 
singular Coulomb term. Further, we shall limit ourselves to case of the strong coupling regime, i.e. when It 
is a well known fact that, in this limit, one can neglect the quantum fluctuations of the phonon field (see p^ j and a 
recent discussion in |p2[). After introducing the canonical transformations 

Ck = + u k , c£ = b£ + it£ (6) 
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and neglecting the quantum components, we then have 

L = ~\rA + 2*" (™) V2 I (re^u k + re"^<) + r <u k - (7) 

(here, all measurements are in the system of units whereby m = h = lj = 1). The classical components of the field 
can be defined from the condition 

9E dE n 

From equations (5), (7) and (8) we obtain 

Here (...) represents the average (*|...|*), where \\&) is the polaron eigenfunction. Substituting (9) into (7) and than 
integrating over k we have 



(10) 



where 

x '\ = V* s (<r\ ) st Vt , x'l = C s (v\) s t(t- 

The emergence of the terms (r) and (r) 2 is conditioned by the changes in the condition of normalization of wavefunc- 
tions in the £-space. Further, we shall consider only s-states of the polaron. The general case will be treated at the 
end of this paper. Given the fact that the eigenfunctions of s-states spherically symmetric, i.e. they depend only on 
r, we can write VE'(^) = vE'(r). Averaging (10) over the angles and using the formula 



i i 



d-r- = — [r'0(r - r 1 ) + r0(r' - r)l, 

^ r 2 + r /2 _ 2x rr> rr /l V ' V > 1 ' 

where the Heaviside function 9(x) vanishes for x < 0, is 1/2 for x — and unity for x > 0, we have 



-\ rA-V2a-^- [ d 4 r] \V(n)\ 2 [r'6(r - r') + r6(r' - r)}+ 
2 \r) J 



+ VzW r I d " v I !^C^7)I :2 I 1 ^CC)I 2 [^'^C^" — ^0 -•- ^"^C^-' — ^'0] — ^-^7^) i^c^)> = o. (ii) 

It is clear that equation (11) is an integro-differential equation of the Hartree-Fock type. This equation can be 
simplified using new units of energy and length, such units are equal to the old ones multiplied respectively by a 2 
and I /a. In this units the parameter a disappears in the last equation. Therefore, we can further put a = 1. 
Let assume that the polaron eigenfunction has the form 

|*) = ^C» (12) 

n 

where 

I / \ n 

yn!(n + 1)! v 7 

It is clear that we need to calculate the matrix elements of the operators constituting L, such as 

(■mi,m 2 \[r 2 0(r 1 - r 2 ) + r\6{r 2 - r 1 )]\n 1 ,n 2 ) = W mim2;ni „ 2 . (13) 
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To this effect we use the integrai representation 

„ +00 +00 

r 2 6{r l -r 2 )+r 1 6{r 2 -r 1 ) = f d Ql f dq 2 e i9iri+i9ar V(gi , q 2 ) , (14) 



2- 



-OO —OO 



where 

00 00 

<p(qi,&)= I drt /dr 2 [r 2 % 1 -r 2 )+r 1 0(r 2 -r 1 )] e -^ ri -^ r2 . 







The operator 



exp(igr) = exp ( ^(2 + N + M + M + ] 



2uj 

can be expressed in the normal form as follows (see, for example, fl9fl , 

exp(igr) = - — - — — cxp( M + J exp( — iVln(l — /i) J exp( M 



(1 — /i) 2 \1 — fi J \ J \1 — fj, 

where, fi = As a result, the calculation of matrix elements of the operator exp(iqr) using the algebra (4) does 
not pose much difficulties. Here, we give only the result obtained 



nil " I ..\m+n— 2s 

(m|exp(^)|n) = ^ -±- g C » n C^ ^ u)m+n+2 , (15) 

where 



(1 - u)™+»+2 ' 



C v „ 



q pKq-pV- 

Subtituting (15) into (13) and (14), after integration over the variables q\ , q 2 , 1*1 , r 2 we finally have the formula 



W _ 1 ( 1 \m l +m 2 +tn+n 2 / (ni + 1) (n 2 + 1) y^^, Sl<oSl + l 

mim2;ni " 2 ~ 16a/ j V (mi + l)(ma + 1) ^ 



si=0 



mi+rai — 2si ti 2 m2+M 2 — 2s 2 

1\*2/^*2 

-2s 2 



X 



ti=0 s 2 =0 t 2 =0 

x (2«i + ti + 2s 2 + t 2 + 3)(2ai +h+ 2s 2 + t 2 + 4)C 2 2 | i 1 + t i 1 + 2 1 S2+t2+2 

'2«i+*i + l / 1 \U/TU 2s 2 +*2 + l / ) 

V ^ °2 Sl +ti + l \p I °2s 2 +t2 + l I (10) 

^ (2s 2 +t 2 +3 + w)2 2s 2+*2+« (2 Sl +tj +3 + u)2 2s i+*i+ M [ ' 1 j 



Other matrix elements are calculated in the same way. 

Let us now consider the zeroth approximation of the OM. The free parameter to is chosen in such a way that the 
condition 

dEW _ 
du 

is satisfied || . Whence 



(0) _ Wnnnn ^(0) _ 



2\/2(n + l) 2 ' " 16(n + l) 4 
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For some first s-states these equations give the following results 

= -0.097656a 2 , e[ Q) = -0.0226 a 2 , E^ = -0.0099 a 2 , 

which are different from the exact ones approximately by 10%. 

Equation (11) is a special case of the following generalized equation (see refs. [9,12]) 

A-E n B^ n )=0. (17) 

Therefore, we describe here a general iteration scheme for the last equation. We shall find the eigenfunction in the 
form 

\^n) = \n) + J2 C kn\k). (18) 

The substitution of this expression into equation (17) gives the system of equations (s is a number of iterations) 

E { n s) = (s»„ + ]T c^B nk ) 1 (A nn + ciUnk) , (19) 

c% = -(A ii -Ei*-VB^~ 1 £ {A ik -E^- 1] B ik )ct l \ (20) 

where the eigenvalue E n and the coefficients of the eigenf unctions Ckn are defined by the relations 

E n = lim c kn = lim c^. (21) 



Equations (18)-(21) completely define quite a simple algorithm for solving equation (17). Apparently, the non- 
diagonal matrix elements of the Hamiltonian arc "included" in each step of the subsequent iteration. It should be 
noted that there exist some other iteration schemes, for example, an interesting variational-iterative one, proposed by 
Burrows and Core (for more details, see (24j|). Here we use the above scheme (18)-(21) because of its simplicity and 
transparency. Moreover, this and some other similar schemes proved to be well convergent in many applications of 
the OM (see [8,9,12] and references therein). 

The results of our calculation show good convergence of the OM in quite a wide range of variation of the parameter 
lu. For s = 18 we obtain the following value for the polaron ground-state energy 

E a = -0.1085128052 a 2 . 

This result is in good agreement with the numerical solution, obtained by Miyake [ p5[ and confirmed by Adamowski 
et al [^6| . It should be noted that good numerical results have been obtained recently by a variational method, based 
on the coherent state representations [^] . We would like to point out that the OM enables one also to find with high 
accuracy the energy and eigenfunction of the excited states, which usually are obtained with much more difficulty 
using other methods J2f|- pij| . We present here the value of the energy for 2s- and 3s-states (for s = 18) 

Ei = -0.02053101 a 2 , E 2 = -0.0083506 a 2 . 

Our results are displayed in table 1 together with the known results for a polaron in the strong coupling limit. 

Table 1. Polaron energy (in a 2 TiujQ units) for s-states in the strong coupling limit. 



Authors 




Is 


2s 


3s 


Pekar (by Miyake 


|5j) -0.108504 






Miyake |2j 




-0.108513 






Adamowski et 


il [ 


|§ -0.1085128 






Feranchuk et al G 


-0.1078 






Smondyrev [p2j 




-0.109206 






Efimov et a^llf 




-0.10843 






Ganbold et al [ 


16| 


-0.107766 






Chen et al |27[] 




-0.10851 






Hagen et al []30| 






-0.02048 


-0.00804 


Balabaev et al 


31 




-0.0206 


-0.00832 


Our results 




-0.1085128052 


-0.02053101 


-0.0083506 
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III. THE OPERATOR METHOD IN SOLVING THE HARTREE-FOCK EQUATION FOR HELIUM 



Let us now consider the problem of helium. For the ground state, we have the following Hartree-Fock equation 
(see, for example, ref. p3|) 



2 n 



ki - r'l 



7 dr' 



*(n) = 



(for the second electron, the corresponding equation is written in analogy with the first). 
The energy of the atom is 



E = 2Ei 



* 2 (n)vE< 2 (r 2 ) 



|ri-r 2 | 

The substitution of the last expression into equation (22) gives 



dY\dY2- 



2 r 



* V) j i 1 f f * 2 (ri)* 2 (r 2 ) 



|r-r' 2 



E 

-dT 1 dr 2 - — 
ri-r 2 2 



*(r) = 0. 



In £-space the equation (23) has the form 



J --rA - Z + -L l dS ^ 2 ( V )\r'9(r - r') 
[ 2 (r) J 



+ r0(r' - r)}- 



(22) 



(23) 



~1(W I / ^ * 2 W* 2 (0[r'd(r" - r') + r"6(r> - r")] + E^j = 0. (24) 

The last equation has the same form as equation (11), for which we have already calculated the matrix elements. The 
key difference is that now we have a negative "coupling constant" a. In order to have the needed solutions, we need 
only to change the formulae of the iteration scheme. Consequently, for s = 18 we have (in atomic units) 

E HF = -2.861679995 

which should be compared with the known result |33| E HF = —2.863 and with a recent calculation |34|] E HF = 
—2.861680. It follows that the OM gives the solutions with very high accuracy independently of the repulsive or 
attractive nature of the interaction. 



IV. CONCLUSIONS 



In this paper we have applied an algebraic approach, namely, the OM to solve the non-linear equations, which appear 
in the problem of polaron in the strong coupling limit and the helium atom in the Hartree-Fock approximation. The 
resuls are as good as the best ones in the literature. Furthermore, we would like to point out that the present approach 
should be successfully applied to the other polaron problems, such as bound polarons and polarons in external fields 
in 3D-space as well as in the other dimensions. 

Finally, we note the following two circumstances: 

(i) As shown in ref. fl35| , the equation, which describes the motion of an individual particle in a system of large 
number of gravitational bosons in fact coincides with equation (11). Therefore, the solutions obtained here for the 
polaron are simultaneously solutions for the above-mentioned problem. 

(ii) We need not here only consider s-states. To this effect it is conducive to use in equation (10) the following 
integral representation 

1 1 f 27r 1 

Further calculation of matrix elements using this formula and the algebra so(4, 2) of operators (3) does not pose much 
difficulty and is affected in analogy with the calculations carried out above. 



6 



[1] Feranchuk I D and Komarov L I 1982 Phys. Lett. A 88 212 

[2] Feranchuk I D and Komarov L I 1982 J. Phys. C 15 1965 

[3] Feranchuk I D and Komarov L I 1984 J. Phys. A 17 3111 

[4] Gerry C C and Silverman S 1983 Phys. Lett. A 95 481 

[5] Witchel W 1983 Phys. Lett. A 97 315 

[6] Feranchuk I D, Fisher S I and Komarov L I 1984 J. Phys. C 17 4309 

[7] Feranchuk I D, Fisher S I and Komarov L I 1985 J. Phys. C 18 5083 

[8] Fernandez F M and Castro E A 1985 Phys. Lett. A 111 104 

[9] Fernandez F M, Castro E A and Meson A M 1986 Mol. Phys. 58 365 

[10] An C Z, Feranchuk I D and Komarov L I 1987 Phys. Lett. A 125 123 

[11] Feranchuk I D, Kuz'min V S and Ulyanenkov A P 1991 Chem. Phys. 157 61 

[12] Le Anh Thu, Hoang Quoc Khanh and Komarov L I 1994 Int. Workshop "Quantum Systems: New Trends and Methods", 

Minsk, May 23-29, unpublished 

[13] Emin D 1989 Phys. Rev. Lett. 62 1544; 1992 Phys. Rev. B 45 5525 

[14] Smondyrev M A, Devreese J T and Peeters F M 1995 Phys. Rev. B 51 15008 

[15] Efimov G V and Ganbold G 1991 Phys. Stat. Sol. B 168 165 

[16] Ganbold G and Efimov G V 1994 Phys. Rev. B 50 3733 

[17] Mitra T K, Chatterjee A and Mukhopadhyay S 1987 Phys. Rep. 153 91 

[18] Lakhno V D and Chuev G N 1995 Sov. Phys. Uspekhi 38 273 

[19] Komarov L I and Romanova T S 1982 Izv. Akad. Nauk BSSR (ser. Fiz.-Mat. Nauk) 2 98 (in Russian); 1985 J. Phys. B: 

At. Mol. Phys. 18 859 

[20] Le Van Hoang and Viloria T 1992 Phys. Lett. A 171 23 

[21] Le Van Hoang, Ly Xuan Hai, Komarov L I and Romanova T S 1992 J. Phys. A: Math. Gen. 25 6461 

[22] Lieb E H and Thomas L E 1997 Commun. Math. Phys. 183 511; ibid. 188 499 

[23] Wilcox R M 1967 J. Math. Phys. 8 962 

[24] Burrows B L and Core P W 1984 J. Phys. A: Math. Gen. 17 559 

[25] Miyake S J 1975 J. Phys. Soc. Jap. 38 181 

[26] Adamowski J, Gerlach B and Leschke H 1980 Phys. Lett. A 79 249 

[27] Chen Q, Fang M, Zhang Q, Wang K and Wan S 1996 J. Phys.-.Condens. Matter 8 7139 

[28] Gabdoulline R R 1993 Dok. Ak. Nauk RAN 333 23 (in Russian) 

[29] Chuev G N and Lakhno V D 1993 Phys. Rev. B 48 12655 

[30] Hagen D E, Van Zandt L L and Prohofsky E W 1970 Phys. Rev. B 2 553 

[31] Balabaev N K and Lakhno V D 1980 Sov. Theor. Math. Phys. 45 936 

[32] Smondyrev M A 1991 Physica A 171 191 

[33] Bethe H A and Salpeter E E 1977 Quantum Mechanics of One- and Two Electron Atoms (New York: Plenum) 

[34] Styszynski J and Karwowski J 1988 J. Phys. B: At. Mol. Opt. Phys. 21 2389; Fischer C F 1987 Comput. Phys. Commun. 
43 335 

[35] Komarov L I 1984 Izv. Akad. Nauk BSSR (ser. Fiz.-Mat. Nauk) 3 90 (in Russian) 



7 



